Approximate locality for quantum systems on graphs 
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In this Letter we make progress on a longstanding open problem of Aaronson and Ambainis [Theory 
of Computing 1, 47 (2005)]: we show that if A is the adjacency matrix of a sufficiently sparse low- 
dimensional graph then the unitary operator e %tA can be approximated by a unitary operator U(t) whose 
sparsity pattern is exactly that of a low-dimensional graph which gets more dense as \t\ increases. Sec- 
ondly, we show that if U is a sparse unitary operator with a gap A in its spectrum, then there exists an 
approximate logarithm H of U which is also sparse. The sparsity pattern of H gets more dense as 1 /A 
increases. These two results can be interpreted as a way to convert between local continuous-time and 
local discrete-time processes. As an example we show that the discrete-time coined quantum walk can 
be realised as an approximately local continuous-time quantum walk. Finally, we use our construction to 
provide a definition for a fractional quantum fourier transform. 



In physics the word locality admits many possible inter- 
pretations. In quantum field theory and condensed matter 
physics locality is understood as the clustering of correla- 
tions iffl H H fi HI III GJ] ■ In quantum information theory 
the quantum circuit model 0, [U reigns supreme as the fi- 
nal arbiter of locality where it is natural to define the non- 
locality of a physical process to be the minimal number of 
fundamental two-qubit quantum gates required to simulate 
the process up to some prespecified error [19]. The cen- 
tral role the quantum circuit model plays in assessing the 
nonlocal "cost" of a physical process strongly motivates us 
to quantify the relationship between the notions of locality 
accepted in other branches of physics and the quantum gate 
cost of quantum information theory. 

Thus, we appear to have at least two different interpre- 
tations of the word locality for quantum systems: on one 
hand we have the clustering of many-particle physics, and 
on the other we have the gate cost of quantum information 
theory. From a physical perspective it is in intuitively clear 
that there should be a strong relationship between these two 
definitions. After all, dynamical clustering implies abound 
on the speed of information transmission. Indeed, for many 
particle systems there are now results quantifying this re- 
lationship: a low-dimensional system which exhibits dy- 
namical clustering can be simulated by a constant-depth 
quantum circuit 19fl. 

However, for the case of scalar and spinor quantum parti- 
cles hopping on finite graphs an explication of the connec- 
tions between the clustering-type interpretation of locality 
and the quantum-circuit type interpretation has yet to be 
completed. An investigation of the graph setting was initi- 
ated by Aaronson and Ambainis lllOll . who established the 
canonical analogues of the clustering and quantum-circuit 
definitions of locality. The main questions remaining are 
now to quantify the relationship between the quantum- 
circuit locality (what Aaronson and Ambainis call "Z- 
locality" and "C-locality") and the clustering-type inter- 
pretation (called "i7-locality") for these systems. 

The objective of this Letter is to provide a (not neces- 
sarily optimal) equivalence between the notions of local- 



ity introduced in El, thus partially resolving one of their 
longstanding conjectures: we establish that a graph-local 
continuous-time process can be written ("discretised") as 
a product of discrete-time processes ("quantum gates"). 
Conversely, we show how to compute an approximately 
local logarithm for a unitary gate which is local on some 
graph G. In other words, we show how to construct a local 
continuous-time process C associated with a local discrete- 
time quantum process T> such that T> may be realised by 
sampling C at appropriate intervals. 

All the quantum systems we consider in this Letter are 
naturally associated with a finite graph G = (V, E), where 
V is a set of n vertices and E a set of edges. We write 
v ~ w if (v, w) E E. We summarise this connectivity 
information using the adjacency matrix A, which has ma- 
trix elements given by A V:W = 1 if v ~ w and A ViW = 
otherwise. For two vertices v,w £ V we let dist(t>, w) 
denote the graph-theoretical distance — the length of the 
shortest path connecting v and w, with respect to the edge 
set E. Let M G M n (C) be an n x n matrix. The spar- 
sity pattern A of M is the n x n {0, l}-matrix given by: 



A 



3,k 



if Mf )fe = and A 



hk = 1 if M hk ^ 0. It 
is sometimes convenient in the sequel to arbitrarily assign 
directions (arrows) to the edges of G. In this case we write 
e + (respectively, e~) for the vertex at the beginning (re- 
spectively, end) of e. Finally, we denote by D(G) the max- 
imum degree of G, which is the maximum number of edges 
which are incident to any vertex in G. 

There is a canonical way to associate a Hilbert space 7i v 
with a finite graph G with vertex set V: we use vertices to 
label a basis of quantum states, so that 7i v = (\v) \ v £ V) 
— this is the Hilbert space of a scalar quantum particle 
constrained to live on the vertices of G. 

We now recall the definitions of locality introduced by 
Aaronson and Ambainis [10] for a quantum particle on a 
graph. Note that the definitions we present here are not 
as general as those introduced in lioll : Aaronson and Am- 
bainis include the possibility of an extra internal degree of 
freedom. While, for clarity, we ignore this extra internal 
degree of freedom it is straightforward to extend our re- 
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FIG. 1: Smoothing out the sawtooth to improve convergence of fourier coefficients: in the first figure the original sawtooth wave is 
shown in red, along with the smoothed version in blue. In the second figure the absolute values of (2iv times) the fourier coefficients of 
the original sawtooth are shown. Note the slow convergence of this sequence. Finally, in the third figure the fourier coefficients of the 
smoothed sawtooth are shown. Note the improved convergence of the sequence. 



suits to cover the more general case. 

Definition 1. A unitary matrix U is said to be Z-local on 
G if Uj.k = whenever j ^ k and (j, k) E. 

Definition 2. A unitary matrix U is said to be C-local on 
Gif: 

1. the basis states \v) can be partitioned into subsets 
Vi, V2, ■ ■ ■ , V q such that Uj^ = whenever \j) and 
\k) belong to distinct subsets P ; ; and 

2. for each j, all basis states in Vj are either from the 
same vertex or from two adjacent vertices. 

Definition 3. A unitary matrix U is said to be i/-local on 
G if U = e lH for some hermitian matrix H with ||-H~||oo < 
7r such that Hj k = whenever j 7^ k and (j, k) E. 

The first result we prove in this Letter shows how an 
H-\ocaX unitary operator may be written as a product of 
C-local unitary operators. This result is entirely standard 
and is a straightforward corollary of the sparse hamiltonian 
lemma of lllll . We sketch a proof for completeness. 

Proposition 4. Let H be the adjacency matrix of a finite 
graph G. Then e ltH may be approximated by a product 
of c\t\D(G) C-local unitary operators, where c is some 
constant. Because a product of C-local unitary operators 
is Z -local on some graph related to G an H -local unitary 
operator is approximately Z-local on some graph G' re- 
lated to G, which gets denser as \t\ increases. 



Proof. The idea behind the proof is as follows. We first 

where hj = Yle^ a * j) 



write H 



V Z5(G) 
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+ h.c, and [ag' , 



a 



(this de- 



composition follows from a colouring of the edges pro- 
vided by Vizing's theorem II 1211 : we denote by Cj the set 
of edges with the same colour). Then we use the Lie- 
Trotter formula to approximate e ltH by powers of Us = 



a iSh„ 



i8h n 



i8h\ s 



\U 



iliii 

L - 2 ;< J 



itH I 



< Q(mA5 + mA 3 \t\5 2 ), (1) 



where m = D(G) + 1 and A = maxj ||/ij||oo — 2, 
where the inequality for \\hj\\ follows straightforwardly 
from, for example, Gersgorin's circle theorem II 1 3fl - Fi- 
nally, we observe that e lShj is a C-local unitary operator, 
for each j = 1, 2, . . . , m. □ 

Proposition 5. Let U be a unitary matrix whose sparsity 
pattern A is the adjacency matrix of a digraph G. If the 
arguments 9 of all of the eigenvalues e lB of U satisfy 6 G 
[0,27r) \ (a, [3), with A = \a — (3\, then there exists a 
unitary matrix V which is H-local on a graph G' given 
by the sparsity pattern of A(G) k where k = c/ (e 2 A), for 
some constant c, such that \\U — V\\oo < e. 

Proof. We begin by writing U in its eigenbasis: 

n 

U = J2^\j)(Jl ( 2 ) 
3=1 

where \j) are the eigenvectors of U and we choose (f)j G 
[0, 271"). By multiplying by an overall unimportant phase 
e^I we can set the zero of angle to arrange for a gap in the 
spectrum spec(C/) of U to lie over the origin. Such a gap 
always exists for finite dimensional unitary operators, but 
not necessarily for infinite operators. 

We want to find a hermitian matrix H so that U = e lH . 
We call this the effective hamiltonian for U. One such 
hamiltonian is simply given by 



ff = X>li>0"l- 



(3) 



While this expression is perfectly well-defined, it is very 
hard to see any kind of sparsity/local structure in H. To 
overcome this we'll find an alternative expression for H 
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FIG. 2: The coined quantum walk [Zona ring with 20 sites. The first figure shows the sparsity pattern of U. The second figure shows 
the spectrum of U with a gap clearly evident. The third figure shows the absolute values of the matrix elements of the logarithm H of 
U constructed in the proof of Proposition [3] Note the rapid decay of the matrix elements away from the original nonzero entries of U. 



defined by Eq. @ as a power series in U. To do this we 
suppose that 



(4) 



and we solve for the coefficients we equate the coeffi- 
cient of \j) (j\ on both sides to find 



4> 3 = c fc e 

k— — oo 

Hence, if we can find such that 



ik(j)j 



o = y] c fc < 



:>k() 



(5) 



(6) 



k=—oo 



for all 9 G [0, 2ir) then we are done. (Recall that we've 
arranged it so there are no eigenvalues of U on the point 
9 = 0.) To solve for c k we integrate both sides of Eq. © 
with respect to 9 over the interval [0, 2tt) against ^e~ lW , 
for I G Z: 

Thus we learn that the c fe are nothing but the fourier coef- 
ficients of the periodic sawtooth function f(9 + 2ttI) = 9, 
9 G [0, 271"), I G Z: 



c fe 



(8) 



'tt, fc = 

,fe> 

Now we know the formula for c k we substitute this into 



Eq. ©: 



H = Yj ( 7r<5 *.o + 



«(1 ~ ^fc,o) 'j fjfe 



(9) 

\ /»• 

fc— — oo 

and truncate the series at some cutoff k < K . If we assume 
the sparsity pattern of U describes a sufficiently sparse 
graph G then U k will also describe a sparse graph for any 



constant k [20], and, as a consequence, the truncated series 
representation for H would also describe a sparse graph. 

Unfortunately we cannot do this: the sawtooth wave has 
a jump discontinuity and hence the fourier series is only 
conditionally convergent. Thus it is impossible to truncate 
the series without a serious error. 

The way to proceed is to assume that we have some fur- 
ther information, namely, that U has a gap A in its spec- 
trum. The eigenvalues of U lie on the unit circle in the 
complex plane so what we mean here is that there is a con- 
tinuous arc in the unit circle which subtends an angle A 
where there are no eigenvalues of U. We arrange, by mul- 
tiplying by an unimportant overall phase, for this gap to be 
centred on the origin. 

The idea now is to exploit the existence of the gap to pro- 
vide a more useful series representation for H. We do this 
by calculating the fourier coefficients of the sawtooth 
wave f(9) convolved with a sufficiently smooth smearing 
function x<y the fourier series then inherits a better con- 
vergence from the smoothness properties of the smearing 
function. That is, we define d k to be the fourier coefficients 
of 



<?(0) = (/*X 7 )(0) 



f(0-y)x-y{y)dy. (10) 



We choose X-yiv) to De a symmetric C°° bump function 
with compact support in the interval [—7, 7] (see the Ap- 
pendix for further details.) Note that, as a consequence 
of the compact support of X-y(y)> = /($)> G 

(7, 2tt — 7). An application of the convolution theorem 
then tells us that the fourier coefficients d k are given by 



(11) 



where ^(co) is the fourier transform of x^iv)- (See Fig.Q] 
for an illustration of the smearing of the sawtooth wave.) 

Using the fourier coefficients d k it is possible to con- 
struct a logarithm J of U which is manifestly sparse if U 
is. We begin by constructing the following approximate 
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hamiltonian: 

oo 

J= S XrW^*- (12) 

fc — — CO 

Choosing 7 < A allows us to conclude that, in fact, H = 
J, because both f(9) and g(9) agree on the spectrum of 
U. 

Our final approximation J k to H is defined by 

k 

Jk = Yl XiUfoU*. (13) 

j=-k 

If U is sparse, with only, say, polynomially many entries 
in n in each row, then so is for j constant. Thus, if we 
choose A: to be a constant, then J k will only be polynomi- 
ally less sparse than U. 

How big do we have to choose fe? To see this we bound 
the difference between H and J k via an application of the 
triangle inequality: 

\\H- JfcHoo < lx 7 (i)l|cil- (14) 

lil>fe 

Now, according to the properties of compactly supported 
C°° bump functions described in the Appendix, XyU) nas 
a characteristic width of I/7, after which it decays faster 
than any polynomial. Thus, choosing k > l/(e J A), for 
an Y j > lj is sufficient to ensure that \\H — Jk\\oo can be 
made smaller than any prespecified accuracy e. 

Now to conclude, we define V = e lJh and use the upper 
bound for 1 1 H — J k \ \ x which we ' ve derived above to bound 
— 

\\U~ V||oo < JfcHoo- (15) 

□ 

Remark 6. By choosing the smearing function x-y (y) to be 
a gaussian a slightly better error scaling can be achieved 
at the expense of a slightly more complicated argument: 
in this case J doesn't equal H and one must bound the 
difference between them. 

Example 7. Consider the coined quantum walk on the 
ring of n vertices: this is the unitary matrix U defined by 
U = (|0)(0|<8>T+|l)(l|<g>T t )iI<g>I,whereTis the unit 
translation operator T\j) = \ j + 1 mod n) and H is the 
hadamard gate ( \ _\ ). The spectrum of U straightfor- 
ward to calculate using a fourier series [ 14]; one finds that 
the eigenvalues A k of U are given by 

Clearly there is a gap A in the spectrum for all n sub- 
tending an angle of 9 with 

9 > 2tan~ 1 (l) = vr/2. (17) 
Thus we find that there exists a logarithm H of U which 
can be expressed as a sum of a few powers of U. Because 
U is sparse, so is H. (See Fig. [2]for an illustration of the 
logarithm of the coined quantum walk.) 



Remark 8. The quantum fourier transform [7, 8] is the uni- 
tary matrix Q defined by the discrete fourier transform: 

n 1 271-ijfc 

Q jtk = -=e » . (18) 

The eigenvalues of Q are well known: because Q 4 = I the 
eigenvalues are the fourth roots of unity. Thus Q possesses 
a gap of size A = tt/2 in its spectrum so we can construct 
a logarithm F of Q as a series Eq. (fT2ll in Q. Although 
F will be dense, it admits a description which is compact 
(i.e., we can efficiently evaluate the matrix elements of F). 
Given the logarithm F it is straightforward to compute the 
square root of Q: y/Q = e^ F = Y^L 4]\ F ' 3 ■ 

Remark 9. Our proof of Proposition [5] also holds for uni- 
tary operators U which are only approximately Z-local, 
i.e., when the condition that C/ j fc = when (j, k) E 
is replaced with Uj, k < e - fidlst (^' £ ) ; or similar. 

The are several questions left open at this point. Per- 
haps most interesting is the question of how to provide a 
combinatorial characterisation of unitary operators which 
possess a gap in their spectrum. Presumably such a char- 
acterisation would take the form of a necessary condition, 
not unlike the isoperimetric inequality II 1 511 - 
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PROPERTIES OF SMOOTH CUTOFF FUNCTIONS 

In this Appendix we briefly review the properties of 
compactly supported C°° cutoff functions. 

Of fundamental utility in our derivations is a class of 
functions known as compactly supported C°° bump func- 
tions. These functions are defined so that their fourier 
transform x 7 (w) is compactly supported on the interval 



[—7, 7], and equal to 1 on the middle third of the interval. 
Such functions satisfy the following derivative bounds 



(19) 



for all j with the implicit constant depending on j. (If we 
have two quantities A and B then we use the notation A < 
B to denote the estimate A < CB for some constant C 
which only depends on unimportant quantities.) This is 
just about the best estimate possible given Taylor's theorem 
with remainder and the constraints that x~/(u>) is equal to 1 
at ijj = and X-y(Lu) is compactly supported. 

The function X-y{t) nas support throughout M but it is 
decaying rapidly. To see this consider 



Xj(t) 



1 

2^ 



-e~^-^-x 7 (u)dLU (20) 
it auj 



which comes from integrating by parts. Continuing is this 
fashion allows us to arrive at 



X-y(u)du (21) 



Since X-y(uj) has all its derivatives bounded, according to 
(fl9t , and using the compact support of Xy(cu) we find 



l*r(t)l < 



it 



- itu *y-1du 



P 1 

Jo MP 
<^ 



(22) 



for all j £ N. Thus we find that x-y (t) decays to faster 
than the inverse of any polynomial in t with characteris- 
tic "width" 1 /j. The existence and construction of such 
functions is discussed, for example, in iflil 171. 



